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We apply the coherency tensor formahsm to the calculation of the spectral distortions imprinted in 
the intensity and polarization of the cosmic microwave background radiation due to the kinematic 
and thermal Sunyaev-Zeldovich effects (SZE). We obtain the first relativistic corrections to the 
intensity produced by the kinematic and thermal SZE, and the first correction to the polarization 
magnitude due to electron thermal motion. 

^ ! PACS numbers: 

o : 

I In addition to the thermal and kinematic SZE which distort the CMB intensity, a CMB polarization signal can 
O . be generated in clusters of galaxies via Compton scattering. The basic process responsible for the generation of 
' polarization is Thomson scattering (low energy Compton scattering) of a radiation field with a quadrupole anisotropy 
' 01 1^ ■ There are several means by which this anisotropy may be generated, in the case of the CMB radiation 
^■f^ ' incident on a galaxy cluster: the primary CMB temperature quadrupole at the cluster, the kinematic quadrupole 
arising from the Doppler boost of the isotropic CMB into the electron rest frame, and double scattering of the 
T-H ■ anisotropic radiation field due to the single scattering thermal and kinematic effects. In this paper we consider only 
^ ' the effect due to the kinematic quadrupole induced by electron motion. 

> Calculation of the frequency dependence of this effect requires a formalism for the treatment of the Compton 
' scattering of a polarized radiation field. In looking at the details of these calculations it becomes apparent that 
, the Stokes parameter formalism conventionally used in polarized radiative transfer and in the primary CMB 
calculations [l^ . is very cumbersome for this purpose, due to the fact that a separate set of polarization basis 
, vectors has to be specified for every photon. Since Compton scattering involves a relativistic scattering electron in 
■ general, Lorentz transformation of the Stokes parameters is necessary, which turns out to be complicated. We were 
[ thus motivated to develop a more elegant formalism for dealing with the Compton scattering of polarized photons, 
O ^ described in 6] (hereafter referred to as Paper I). 

Previous calculations have determined the intensity 0, ^3 1 ^^'^ polarization magnitude 0, S S IS Q| of the 
distortion of the scattered radiation field as expansions in the dimensionless electron temperature 9e = ksTe/meC^ 
and dimensionless bulk velocity magnitude Pb = Vb/c with various formalisms, but not in such a systematic and 
^ , explicit fashion as we describe here. We present a detailed calculation of the polarization matrix of the scattered 
• radiation, which yields in addition to the polarization magnitude, the unpolarized thermal and kinematic effects also. 
. ^ ' In the intensity we obtain the first relativistic correction to the thermal and kinematic SZE, and their cross term. 

In the polarization we obtain the first correction due to thermal electron motion. We work entirely in the Thomson 
?H limit (neglecting electron recoil). 

We break the calculation into two stages. In fjj we perform the calculation in the case of a clump of electrons 
with zero temperature moving with a collective bulk velocity /3f, along the z-direction in the "lab" frame (the CMB 
rest frame), working entirely in the rest frame of the electroirs. The polarization matrix of the scattered radiation is 
obtained, which on transformation to the lab frame yields the kinematic effects to any desired order in In we 
extend this calculation to allow for thermal motion of the electrons. This is done by first generalizing the calculation 
of the rest frame scattered irratrix in f^to the case of a lab frame electron velocity in an arbitrary directioir. Since the 
algebraic manipulations are lengthy and tedious, a computer algebra system is used (one of advantages of our formalism 
is that it is quite simple to implement on a computer algebra system capable of handling matrix manipulations). After 
transformation of the resulting scattered beam into lab, the integration over electron velocities is performed. The 
electroirs are assumed to have a phase space density ge given by a relativistic Maxwellian distribution with electron 
temperature Tg and a bulk 3-velocity VJ, = /3bC. 
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I. COLD ELECTRONS 

In this section we use the formahsm developed in Paper I to compute the CMB intensity and polarization distortion, 
in the approximation of a single scattering, due to scattering of the unpolarized isotropic part of the incident CMB 
intensity from electrons moving with a given bulk 3-velocity Vj, = /3f,c with respect to the CMB rest frame. The CMB 
rest frame will henceforth be called the "lab frame" in this section. We deal only with an idealized galaxy cluster 
composed of a concentrated clump of electrons of density rig and corresponding optical depth tx in lab frame. 

In the single scattering limit, since the resulting scattered radiation field must be symmetric under rotations about 
the electron bulk velocity, the angular dependence of the intensity and polarization magnitude of the scattered 
radiation is a function only of the angle cosine fj, — n ■ f3b/\f3b\ between the electron bulk velocity and the line of 
sight. By symmetry, the actual polarization vectors on the sky produced by this effect are simply all orthogonal to the 
direction of the cluster bulk velocity (they could also be parallel to it, depending on the sign of the Stokes parameter 
Q, but turn out to be orthogonal '9'|). Note that in a calculation with a real cluster with spatially extended structure 
we may replace tt with the optical depth integrated along the line of sight to obtain the intensity distortion for each 
viewing angle. 

For simplicity we choose, without loss of generality, to align Vf, with the z-axis of a Cartesian coordinate system. 
Our task is to calculate the polarization matrix resulting from Thomson scattering of the incident unpolarized CMB 
blackbody radiation in the lab frame into the lab frame viewing direction n. We begin by writing down the intensity 
and polarization matrix of the incident photons in both lab and rest frames. Primes denote the rest frame, unprimed 
quantities denote lab frame. We align the velocity 3-vector of the electrons in lab frame with the z-axis, and write 
the electron velocity in lab frame coordinates as 

Ve = 76(1, PbZ) , 7& = . ^„ . (1) 

V 1 - 0b 

The lab frame 4- velocity is denoted vi. The rest frame momentum of the incident photon is = p'^{l,n[), where 
the rest frame direction vector is expressed in polar coordinates with respect to the z-axis: 



^cosV'-^/l -Aif , sinV'-\/l-/if , 



(2) 



The coordinate system is illustrated in Fig.Q The corresponding lab frame momentum is = pi(l,ni), where the 
lab frame direction vector is: 



rii 



(^COStpiy/l - ^i'^, Sini/'i\/l - Mi^; Mi) • (3) 



Assuming unpolarized isotropic incident CMB radiation in lab frame, the intensity polarization matrix of a photon 
incident in the lab frame with 4-momentum pi is given by 

I^''ip.,vi)^hoipdP''''ip.,vi) . (4) 

Here we restrict attention to the case where Iq is the Planck fimction at the mean temperature of the CMB, Tcmb: 

2 c p^ 

h{Pi) = T^gp.c/fc^TcMB _ 1 ■ 

The incident photon momentum in the lab frame is Doppler shifted on going to the rest frame: 

P't = 7fcP»(l - Pb^J■^) , (6) 
This may be written in terms of the incident polar angle in the rest frame. Using the formula for relativistic aberration, 

1 + PfcMi 



we obtain 



76(1 + Pb^J-'^) 



(8) 
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FIG. 1: The coordinate system in the (a) lab frame, and (b) rest frame, used to evaluate the polarization matrix. In lab, the 
clump of electrons indicated at the origin travels along the a-axis with velocity Vb- Note that in the rest frame, we choose to 
consider the photons scattered in a direction n' in the x-z plane, but the incident photon direction n'i is in a general direction. 



The specific intensity tensor in the rest frame can be obtained from that in the lab frame using the transformation 
law of the intensity between frames, 



(9) 



c2 g76Pj(l+/3i'A'i)/'=B^CMB _ \ 



The isotropic specific intensity in the lab frame, Ia{pi), transforms into an anisotropic intensity in the rest frame 
which is still of blackbody form but with a temperature with angular dependence: 

The incident radiation field in the rest frame is of course also unpolarized and has intensity normalized polarization 
tensor: 

In Paper I, the following kinetic equation for the time evolution of the polarization matrix due to Thomson scattering 
in the electron rest frame was derived (Eqn. (168), Paper I): 



^I^"'{p,Ve) = neCTT 

at 



P^iP^ ^e)P%{p, Ve)I''^\p^,Ve) " /^'^ (p, V. 



(12) 



where 



1 T) TJ 

Pi_,„{p,v) = rii_,^ + -{p^v^ + p^v^) ^ where p=-v^^p^. (13) 

p p'^ 

Here we have changed the polarization matrix normalization from phase space density to intensity, which is valid in 
the Thomson limit since the incident and scattered photons have the same energy in the rest frame. The Thomson 
limit form is appropriate here since even with the boost from the lab to rest frame, the CMB photon momenta are 
a tiny fraction of the electron mass and therefore recoil is negligible. We evaluate the kinetic equation at photon 
4- momentum p, with the following components in rest frame coordinates: 



p^' -p'(l,n'), n'- (^cosV;'yi~^,sin^'yW^,Ai'^ , 



(14) 



where p' — p[ since we are working in the Thomson limit. 

In the single scattering limit, we may insert the polarization tensor of the incident unpolarized radiation field in 
the right hand side of the kinetic equation to obtain the scattered beam: 



(15) 
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where 



(16) 



The 00 and Oi components of this tensor equation obviously vanish when evaluated in electron rest frame coordinates. 
We evaluate the gain term by first performing the integral over azimuthal angles dip'^, using the explicit form of the 
spatial part of the projection tensor: 



( 1 + - 1) cos2 V^,; (K'-l)cosV'^sinV'^ -^^y^ l - cos ^ 



(K' - 1) cos sin 1 + (^f - 1) sin2 



Performing the integral over azimuth yields 



/2 



(17) 



^(1 + Mf) 



I pi J 



2tt 








(18) 



i(i-Mf) 



To further simplify, we may evaluate the rest of the scattering term at ■0' = 0, since by azimuthal symmetry 
of the radiation field in the rest frame the intensity tensor for a general {^',tp') is related to that at ijj' = by 
a rotation about the z-axis through angle tp', and the polarization magnitude is independent of azimuth. Putting 

n' — 1 — /z'^, 0, ^' j , the explicit form of the spatial part of the projection tensor is (p, Vg) is: 



-mVi - / 



(19) 



-n'^/l-fi'^ l-/i'^ 



Now performing the multiplication of two of the matrices in Eqn. H19() with the matrix in Eqn. (jlSfl . yields the 
azimuthal integral of Eqn. H16|l required in the gain term of the kinetic equation H15(l : 



-/iVl-/i''G||(/i',/iO 

G±(/i',M:) 



(20) 



where 



- 2 



One can check that 



(21) 



(22) 



(evaluated at V'' = 0). Thus as Pb 0, {p,Ve) 0, since by symmetry scattering of an isotropic radiation field 

from a stationary electron cannot alter the radiation field. 

Now putting together the gain and loss terms, and integrating /i^, we find for a finite rest frame time interval At' 
(all evaluated at ■0' = in rest frame coordinates) : 



AI'{p,v, 



-^l'V^^^{p',^i') 

I±ip',fi') 
-/iVl - (p', m') /|i (p', ^') (i - A^'') 







(23) 
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where we set the optical depth in the rest frame to t!^ = n'^aT^t' , and defined 



(2 - /i") J(P') + (3/^'' - 2)K{p') - ll'ip', 



J{p') + K{p')~^I\p',ti') 



(24) 



and the functions (note that these are functions of p^, but in the Thomson hmit this equals the scattered photon 
momentum p'): 



•Hp) 

K{p') 



-1 
1 



(25) 



Using the dimensionless frequency x' = cp' /ksTcuB, and the constant io — 2{kBTcMB)^ / (hc)'^ , these functions are 
given by the integrals 



p762;'(1+/36mD _ I 



J{x',Pb) = iox'^ J dfi'i 

K{x',f3b) = iox'^ j\'i'dfi'^[e-'^^'^^+''^^'^^ -1 

The intensity of the scattered radiation in the lab frame is thus given by the trace 
A/'(x',m') = TT[AI'ip,v,)] 



(26) 



3 , 



(3 - ^l'')J{x\ Pb) + (3m'' - l)K{x\ M - fl'ix', m') 



where 



I'ix',^i')^iox'• 



,7i,x'(l+/3i,M') _ 



(27) 



(28) 



The polarization magnitude, in the limit of small is given by the formula for the polarization of a perturbed 
unpolarized beam, given in Eqn. (52) of Paper I: 



U{x',f,')^UiI'ip,v,) + AI'{p,v,))^(^^^'^ n(AJ') , 



where 



Evaluating this, we find 



n^(A/') 



2Tr[Ar2] 



Tr[A7']^ 



, _ 3 , , ,2N 5K{x',Pb) - J{x',f3b) 



(29) 



(30) 



(31) 



This formula, and the /i,^ angular dependence inside the integrands of K and J, shows that only the quadrupole in 
the incident intensity in the rest frame generates polarization. The polarization magnitude has the familiar sin^ 9' 
dependence. As (3b 0, the incident radiation field in the rest frame becomes isotropic, I'{x',fi') — + Iq (of Eqn. Q), 
yielding J 2/o, K 2/o/3, and thus n(a;', /i') ^ as expected. 

We now expand the integrands of J{x' , (3b) and K(x' , (3b) in powers of 



no{x) 



(e^' - 1)- 



E 

n=0 



(-1)"5 



n en+1 



(32) 
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where 

J ^ g(7.-l)x'g7->/3.M':i:' _ 1 _ (33) 

Expanding S up to second order in /3b, we find: 



d^l' 5 = I3lx' ( 1 + ^a;' ) + 0(/3^) 



dti' /i'^ 5 = \plx' (l + + 0{l3t) , 



-1 

\^V"<5' = |/3,V2+0(/3,^). (34) 
-1 5 

To obtain Il(x' , to second order in /3t,, since the numerator is already second order, we may replace I'(x',fi') in 
Eqn. lini) with Iq. 

Using the results in the previous equation we obtain finally the polarization magnitude in the rest frame to second 
order in /3t: 

n(x',M') = 2^ ^T^'' 'l}f_^^y (1 - m'^) ■ (35) 

On Lorentz transforming into the lab frame, the polarization magnitude of the scattered photons does not change, 
but the photon angle is aberrated, with 

Since /i'^ — ^"^ + 0{Pb), and x' = x + 0{f3b), and the optical depth rfj, transforms like dt', we have the same result in 
the lab frame quantities to this order: 

n(-,A^) = ^-T.^^^/3?(l-M^)- (37) 



This result was obtained before by 0,0,01 using different methods. Note that the (1 — /i^) dependence implies 
that this component of the CMB polarization is a direct measure of the peculiar velocity of the cluster gas perpendicular 
to the line of sight, which in conjunction with the intensity measurement allows, in principle, measurement of all 
the components of the cluster peculiar velocity. This will be an important cosmological probe, if the polarization 
measurements can be made. This will be a considerable experimental challenge, since the polarization magnitude 
is rather small, typically 0.1 /iK at most, as illustrated in Fig. |3 (note that cluster bulk velocities rarely exceed 
1000 kms~^). In panel (b) we show the polarization magnitude as a Rayleigh-Jeans (RJ) brightness temperature, 
given by 



ATb = V^Ql + AC/2 . (38) 

One might worry that the dimensionless polarization magnitude 11 goes quadratically in a; as a; ^ cxi, which would 
seem to be a problem since the polarization magnitude must be bounded by unity. However, the analysis we have 
given is only the lowest order result - at high photon frequencies, relativistic corrections will modify Eqn. (|37|l . Since 
we essentially expanded in powers of /JfcX, our analysis cannot be trusted for frequencies greater than x w l/(3b- 

We finish this section by expanding the total intensity of the scattered radiation, given in Eqn. (|27|l . to second order 
in (3b, and performing the transformation to lab frame to obtain the kinematic SZ distortion and its first relativistic 
correction. To do this calculation it is convenient to work with the phase space density rather than the intensity. The 
intensity distortion AJ'(x',/i') is related to the phase space density distortion by 

Arix',n') ^ Af'{x',f,') 

/'(x',M') r{x\ti') ■ ^^^^ 
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0.1 1 10 0.1 1 10 



x=hu/K T x=hu/K T 

(a) Dimensionless II/tt (b) Brightness temperature 11 

FIG. 2: Polarization magnitude of the CMB scattered by a concentrated cloud of cold electrons with a bulk flow velocity Vb 
transverse to the line of sight. In (a), we plot the dimensionless polarization magnitude n(x) divided through by the optical 
depth tt- In (b), we plot the polarization magnitude n(a::) as a (Rayleigh- Jeans) brightness temperature distortion, taking 
tt = 0.01. The upper and lower curves in both (a) and (b) correspond to Vj, — 5000 km s~^ and Vj, = 1000 km s~^ respectively. 



From the transformation law of the left hand side of the kinetic equation, given by Eqn. (204) of Paper I, we find an 
equation for the rate of change of the phase space density in lab: 



1 



df{x, ^) 



/(a;,/i) dt 



(40) 



Expanding the functions K{x' , Ph) and J{x' , fib) in Eqn. (|27|) up to second order in (3h as before, and expressing the 
right hand side in lab frame quantities by making the replacements x' = 7ba;(l — /Sb/i), fi' — {fi — [3b) /{I — /3h/i) and 
rig = ne/76, we find to O{01) the fractional intensity distortion in lab: 



A/(x,/i) 
Io{x) 



tt- 



1 



/^/36 - 



, , + coth(f) ^ 
M + 



(41) 



where the lab frame optical depth is defined by tt = neCTxAt. This is the first relativistic correction to the kinematic 
SZ effect, obtained previously by Note that without the correct "flux factor" 7f,(l — /J^/i) in Eqn. (|40|) . this would 

differ at second order in /3^. The first term is simply the lowest order kinematic SZ distortion, where Vr = —fJ-Pb is the 
bulk velocity projected along the line of sight (which is opposite to the direction of the scattered photon momentum, 
hence the minus sign). 



II. HOT ELECTRONS 



We now extend to the more general case of a Maxwellian distribution of electrons with dimensionless temperature 
6e = kBTe/nie uioviug with a bulk velocity Vb with respect to the CMB rest frame (lab frame). In the single scattering 
limit, the Thomson scattering of isotropic blackbody radiation from a Maxwellian distribution of electrons moving 
with a bulk velocity Vb — (3bC produces a scattered radiation field whose intensity and polarization magnitude are 
azimuthally symmetric about Vb- Our goal is to compute the polarization matrix of the scattered radiation field, as 
an expansion in powers of Vb and 9. This computation will yield, to lowest order, the usual thermal and kinematic 
SZ distortion of the intensity, and the polarization magnitude to lowest order in Vb- Going to higher order yields 
the "interference" terms between the thermal and kinematic effects, in both the intensity and polarization, and the 
relativistic corrections. 



8 



The first task is to determine the lab frame polarization matrix of the scattered beam due to scattering of an 
incident unpolarized isotropic blackbody radiation field in lab by an electron with a general lab frame velocity (3. 
This is not the bulk velocity but rather the velocity of some of the electrons in the thermal distribution, which will 
eventually be integrated over. This part is just a generalization of the calculation performed in section ^ The 
resulting polarization matrix of the lab frame scattered radiation field as a function of electron velocity may then be 
averaged over a distribution of lab frame electron velocities to yield the observed lab frame result. The steps required 
to compute this are described below. The actual calculation, even at lowest order, is quite lengthy, so a computer 
algebra system (Mathematica 12]) was used to perform the calculation. We do not give all the algebra but just 
outline the procedure. Henceforth in this section primed indices refer to components of 4-vectors in the electron rest 
frame, and unprimed indices to components in the lab frame. 

It is convenient to integrate over angles in the electron rest frame, but to express the electron velocity and final 
state photon momentum in lab coordinates throughout (to avoid a cumbersome transformation of rest frame angles 
to lab frame). In lab frame coordinates, the electron 4- velocity is 

< = 7(l,/3), 7^_i==, (5^\(3\. (42) 

The Cartesian coefficients of (3 are denoted In rest frame coordinates, the velocity of the lab frame is of course 

<' =7(1,-/3)- (43) 

The scattered photon momentum in lab frame coordinates is written 

p^=p,(l,n,). (44) 

To simplify the computation, we set up a polar coordinate system with polar axis along the z-direction and evaluate 
the scattered polarization matrix at azimuth ijjs = Q'. 



ris = v/l-M?, 0, ^^s) . (45) 



This is no loss of generality provided we choose the bulk velocity Vf, to lie along the z-direction, in which case the 
polarization matrix for a general ris is related to the one calculated here by a simple rotation about the z-axis. 

The scattered photon momentum in the electron rest frame is found by applying Lorentz transformation matrices 
to obtain = A'' ^(/3) p^, where 



AO - 7 = 1/v/T^, A\=A\ = -lf3' , 

= (7-l)^+^%-- (46) 
Using the notation of §??, we denote the momenta of the incident photons in the lab and rest frames as follows: 

Pi = Piil,ni) , =p^(l,nO , 

n', = (^cosV^^W?' sin^/^^Tw?, m^) ■ (47) 

with Pi = ~vi ■ Pi, p'^ = —Ve ■ Pi- In the lab frame, the scalar occupation number of the incident photons is isotropic 
with a Planck spectrum: 

(Note, do not confuse rii, a direction vector, with Ui, the occupation number!). In the rest frame, n[{pi) = ni{pi), 
but the occupation number of the incident photons is no longer isotropic since photons with different momenta are 
aberrated through different angles. Thus rii becomes a function of and through pi: 

KiPt) = — , , ^ , (49) 

where in terms of /3, 

p,(p>'J=7P-(l + /3-<) ■ (50) 
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The angular dependence of the incident radiation field in the rest frame is obtained by expanding H49|l in powers of 
the velocity components f3i. For the lowest order polarization computation, the expansion must be taken up to at 
least second order in the velocity components. 

Then as in the previous section, the right hand side of the rest frame master equation (|12|l is constructed, and 
the integration over the rest frame angles of the incident beam performed. The resulting matrix is then transformed 
into lab frame by the application of two projection tensors, and the lab frame fractional intensity distortion obtained, 
making sure, as in Eqn. (|40|l . to multiply by the correct flux factor, which now has the form 7(1 — (3 ■ n). 

We thus obtain the lab frame polarization matrix as a function of the lab frame photon direction and the velocity 
components (3. In the lab frame, the integration over electron velocities is performed. To do this we need first to 
construct the distribution function of electron velocities in lab frame. In the "comoving frame", denoted with primes, 
in which the average electron velocity vanishes, the electron phase space distribution function as a function of the 
electron 3- momentum q' is assumed to be a relativistic Maxwellian at dimcnsionlcss temperature O^'. 

g^{q')^goe~^('''yMK (51) 



where E{q') — ^ q'^ + m^, and go is a normalization constant which depends on the total number density of electrons. 
We use a relativistic Maxwellian in order to retain the corrections to the SZ effect in a mildly relativistic plasma. 

With q' = me7'/3', where /3' is the electron 3- velocity in the comoving frame, and 7' = 1/ \J\~ we have (as a 
function of /3 since the distribution is isotropic in the comoving frame) 

5e(/?')=5oe-^'/^= . (52) 

The number density of electrons in each comoving frame momentum element Sp' is thus 

d< = 47rgop'^ dp'e-^'/'^'' . (53) 

Integrating this distribution over the element (Pp' yields the electron number density in the comoving frame n'^. Using 
p' = rUeiP', we find p'^dp' = ml^f'^ (3'^ d(3' . Thus 



.90 



4^7'5/3'2d/3'e-T'/^= 



= ATrgomi9eK2il/0e) , (54) 

where K2 is a modified Bessel function (for a derivation of this result see for example ^3). 

Thus in the comoving frame, the number density of electrons in each comoving frame velocity element d^v' , is 

dn' = n' — —— — ■ — - . (55) 

For small 9e, the denominator can be expanded: 

AiT9,K2il/9,)e'^'>^ - (2^0e)'/' [1 + 150e/8 + • • • ] , (56) 

yielding the familiar prefactor of the non-relativistic Maxwellian to lowest order. 

Now we wish to compute the analogous lab frame quantity by a Lorentz transformation from the comoving frame. 
Since the distribution function is a Lorentz scalar, the number density element transforms like the momentum space 
volume element in Eqn. H53f) . Using the Lorentz invariance of the quantity d^p/p = d^p' /p', it follows that 

p'^dp' = mli^(3'^d(3' = ml-i^P^d(3{i /-^) . (57) 

Choosing the bulk velocity of the comoving frame with respect to the lab frame to be Vf, = P^c in the z-direction, we 
have 

ih - 7b (1 - PzM , 7b = (1 - . (58) 

For calculations it is convenient to write the distribution function in lab frame in a form in which the non-relativistic 
part of the Maxwellian, which has Gaussian form, is pulled out and the rest expanded in a series in powers of the 
velocity relative to the dimensionless bulk velocity /3h = (0, 0, /3b): 



dn, e-(^-i)/« 



He 6,e^/<^'K2{l/9, 



iH (1 - Mb) e-^[^^(i-^^'^^)-il/<'el ^ . (59) 

J 477 
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0.1 1 10 0.1 1 10 

x = ht'/K T x = ht'/K T 

(a) Thermal effect (b) Kinematic effect 

FIG. 3: The frequency dependence of the thermal and kinematic SZ effects and their first relativistic corrections (as RJ 
brightness temperature distortions, as defined in Eqn. 13811 '). The solid lines in each plot show the lowest order effect, and the 
short dashed lines show the first relativistic correction. The long dashed line in plot (a) shows the "interference" term Fi^. 



Making the substitution /J^ ^ /3z + VJ,, the part in square brackets may be expanded straightforwardly about unity 
in powers of l3x,Py,Pzi Pb and 9e. Defining 

f3^ = Pl + Pi + M , (60) 

the exponential factor in front can be written as 

g-/3V2e. X prefactor , (61) 

where the prefactor is an expansion about unity in powers of Pz,Px,f3y. The result is a Gaussian multiplied by a 
prefactor which is polynomial in the components /?; with coefficients which are functions of f3b and 6e. 

A further transformation is required before the lab frame integral can be done. The integral (fif3 ranges over the 
velocity sphere |/3| < 1. To simplify the Gaussian integrals, it is easier to make the transformation Pi — > Ui/"/, and 
integrate d^u over all space. With this transformation, we find 

7 = VTT^ , "f^P^dp = u^du . (62) 

Steps similar to those described above yield an expansion about the transformed bulk velocity f3i, ^ Ut = 
(0, 0, Pb/ \/l + Pb) ill powers of Uz = Uz — Pb/ \/i — P^, Ux,Uy. This form is then convenient for integration by a 
symbolic algebra package. 

We expand the prefactor to terms up to sixth order in the coefficients, and up to second order in both Pi, and 
de- Integration over the electron distribution function then yields the lab frame polarization matrix as a function of 
the bulk velocity Pt and electron temperature 6^. Taking the trace of this matrix gave the following result for the 
intensity distortion: 

Iq — 1 

+ PbF^f'i^i) + PlF^{x, fi) + BePhFl'^ix, (!) + ■■■] . (63) 
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FIG. 4: Polarization magnitude generated by scattering of CMB monopole (solid line), and the first relativistic correction 
(dashed line), as RJ brightness temperature distortions, in the case of a concentrated cluster with electrons at temperature 
ksTe = 10 keV, a bulk flow velocity Vb = 1000 km s~^ perpendicular to the line of sight, and an optical depth to scattering of 
TT = 0.01. 



Here Fq is the well known thermal SZ distortion piece, and Ff is the first relativistic correction to the thermal effect: 

Foi^) = -4 + X coth (^1) 
[5(-l + e"^^^ 



{x) 



X (-235 + 77x2) cosh (I) 
3 



+x (235 + 7a;2) cosh 
-8 (-25 + 42x^ + (25 + 21 x'^) cosh(a;)) sinh 



(64) 



The frequency dependence obtained here agrees with that obtained by 

The terms F^^ and F^ are the lowest order kinematic effect and the first relativistic correction respectively: 



F/^(x,M) 



-1-M^ 



x(3 + ll/i|2) coth(f) 
20 



(65) 



which agree with the forms in The "interference" term between the thermal and kinematic effects is: 

[-45 + 14a;2 + (45 + 7 x^) cosh(a:) - 47 a; sinh(x)] 



10 sinh^ (x/2) 



(66) 



The thermal and kinematic effects, their relativistic corrections, and the interference term are plotted for representative 
cluster parameters in Fig. 13 These were computed for a cluster with electrons at temperature ksTe = 10 keV, a bulk 
flow velocity Vb = 1000 km s~^ at an angle cosine = l/v^ to the line of sight, and an optical depth to scattering 
of tt = 0.01. (Note that the dips in the curves are zero crossings). 

Computing the polarization magnitude of the final lab frame matrix, we find: 



where 



u{x,fi) = TTP'b {'^-f^')K{x) + Frix)ee + o{el)] , 

(1 + e^) a;2 



F^{x) 



20(-l + e^)^ 



(67) 



(68) 



12 



and 



e 2 X 



10(-l + e^) 



-4+ llx^) cosh ( - 



8x3 sinh 



(f) 



'3a; 

(4 + x'^) cosh ( — 



'3a; 
sinh ( — 



(69) 



As 0e ^ this reduces to the cold electron result Eqn. (|37|l . The frequency dependence of these results for a cluster 
with typical parameters is shown in Fig. 0] 
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